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Abstract 

Creontiades dilutus (Stal), the green mirid, is a polyphagous herbivorous insect 
endemic to Australia. Although common in the arid interior of Australia and 
found on several native host plants that are spatially and temporally ephemeral, 
green mirids also reach pest levels on several crops in eastern Australia. These 
host-associated dynamics, distributed across a large geographic area, raise ques- 
tions as to whether (1) seasonal fluctuations in population size result in genetic 
bottlenecks and drift, (2) arid and agricultural populations are genetically iso- 
lated, and (3) the use of different host plants results in genetic differentiation. 
We sequenced a mitochondrial COI fragment from individuals collected over 
24 years and screened microsatellite variation from 32 populations across two 
seasons. The predominance of a single COI haplotype and negative Tajima D 
in samples from 2006/2007 fit with a population expansion model. In the older 
collections (1983 and 1993), a different haplotype is most prevalent, consistent 
with successive population contractions and expansions. Microsatellite data 
indicates recent migration between inland sites and coastal crops and admixture 
in several populations. Altogether, the data suggest that long-distance dispersal 
occurs between arid and agricultural regions, and this, together with fluctua- 
tions in population size, leads to temporally dynamic patterns of genetic differ- 
entiation. Host-associated differentiation is evident between mirids sampled 
from plants in the genus Cullen (Fabaceae), the primary host, and alternative 
host plant species growing nearby in arid regions. Our results highlight the 
importance of jointly assessing natural and agricultural environments in under- 
standing the ecology of pest insects. 



Introduction 

Many insects that damage agricultural crops have invaded 
the resources provided by agriculture across wide areas 
and this has generated alternative predictions as to their 
evolutionary trajectories. Although the provision of novel 
resources by agriculture might promote host adaptation 
(Via 1990), it has also been argued that gene flow will 
increase among populations of native insects when their 
range is expanded through the anthropogenic spread of 
potential hosts, making local adaptation less likely (Oliver 



2006). Insects that use native and introduced hosts thus 
provide a "natural experiment" to explore the likely con- 
sequences of ongoing anthropogenic change in plant dis- 
tribution and abundance. 

Few genetic studies have examined the interactions of 
insects between native host plants and agricultural 
resources simultaneously, but the available evidence indi- 
cates that several outcomes are possible, including 
geographic differentiation, host-associated differentiation, 
and widespread gene flow. Both the rice mirid Stenotus 
rubrovittatus (Hemiptera: Miridae), native to Japan, and 
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Queensland fruit fly Bactrocera tryoni (Diptera: Tephriti- 
dae) show strong geographic differentiation (Yu et al. 
2001; Kobayashi et al. 2011). In the former, it indicates 
divergence across Pleistocene refuges and the latter diver- 
gence as invading crop and fruit hosts outside its original 
Queensland distribution. Furthermore, an isolated fruit fly 
population in inland Australia (Alice Springs) showed 
strong genetic evidence of a population bottleneck. Host- 
associated differentiation has also been recorded in the corn 
leafhopper Dalbulus maidis (Hemiptera: Cicadellidae), and 
this has been associated with a shift from wild hosts to 
maize (Zea mays) within the last 9000 years since domesti- 
cation (Medina et al. 2012). In contrast to the above exam- 
ples, a lack of isolation by distance was found across 
1700 km in the migratory moth Trichoplusia ni (Lepidop- 
tera: Noctuidae), between its native range in California, 
and crops that it seasonally invades in Canada (Franklin 
et al. 2010). Clearly, the patterns found to date are strongly 
influenced by the biology and life history of the organism 
in question, as well as the environment it inhabits. 

Strong regional differentiation, as found in Stenotus ru- 
brovittatus (Kobayashi et al. 2011), and Queensland fruit 
fly (Yu et al. 2001), might be expected in species that do 
not regularly disperse long distances (Bohonak 1999). 
Conversely, high gene flow, as documented for the migra- 
tory moth Trichoplusia ni (Franklin et al. 2010), has also 
been reported in many widespread agricultural pests 
(Endersby et al. 2006, 2007; Margaritopoulos et al. 2009), 
and even in pest species thought to be relatively sedentary 
(Voudouris et al. 2012). Anthropogenic and unassisted 
dispersal can both allow the invasion of novel resources 
by insect populations (Stone and Sunnucks 1993; Stone 
et al. 2007). Such anthropogenic dispersal was thought to 
be the primary mechanism allowing colonization of grain 
storages by Tribolium castaneum, as this species was con- 
sidered relatively sedentary (Drury et al. 2009). Active dis- 
persal by flight has subsequently been shown to better 
explain patterns of regional genetic differentiation (Ridley 
et al. 2011; Semeao et al. 2012), highlighting that the 
capacity of organisms to disperse can be underestimated. 

Another important aspect of pest insect dynamics is fluc- 
tuations in population size, which are expected based on 
the seasonal availability of most agricultural crops and the 
occurrence of pest outbreaks. Temporal fluctuations in 
gene frequencies are tied to the number of effective breed- 
ers (Waples and Teel 1990), and the temporal stability of 
regional genetic structure recorded across 5 years of sam- 
pling in the Queensland fruit fly implies that populations 
of sufficient size persist across seasons, despite the occur- 
rence of regional outbreaks of this species (Yu et al. 2001). 
Regional differences in outbreak propensity in the migra- 
tory locust (Locusta tnigratoria) have allowed an empirical 
evaluation of its effects (Chapuis et al. 2008, 2009). No dif- 



ference was found in genetic diversity between outbreak 
and non-outbreak populations (indicating that non-out- 
break populations persist in sufficient size), but regional 
differentiation was much higher for non-outbreak popula- 
tions (Chapuis et al. 2008, 2009). Spatial and temporal var- 
iance in population size, migration rates, and extinction 
rates are predicted to not only affect mean _F S t, but also 
result in large fluctuations in the genetic differentiation 
between populations over time (Whitlock 1992). 

Patterns of host-associated differentiation might be 
obscured by migration, bottlenecks, and population 
expansion, and interpreting the relative effects of demo- 
graphic processes remains a challenge for empirical popu- 
lation genetics (Pavlidis et al. 2008; Li et al. 2012). 
Furthermore, host plant-associated differentiation follow- 
ing host shifts is considered more likely in host specialists 
(Funk et al. 2002). Few studies report host-associated dif- 
ferentiation in insects that use multiple host plants (but 
see Sword et al. 2005). Geographic differentiation gener- 
ally appears to be higher in host plant specialists than 
insects that use multiple hosts (Kelley et al. 2000; Zayed 
et al. 2005; Gaete-Eastman et al. 2007; Habel et al. 2009; 
Groot et al. 2011). This correlation is usually interpreted 
as a consequence of the spatial patchiness of a single 
resource in comparison to the more widespread availabil- 
ity of multiple resources. Resources can, however, be tem- 
porally patchy and it is not clear whether this correlation 
would hold under these circumstances. 

In this study, we examine dispersal, fluctuations in popu- 
lation size, and the use of multiple host plants in Creonti- 
ades dilutus, the green mirid, across its native range (and 
host plants) in arid regions of Australia, and novel hosts 
(agricultural crops) that it has invaded within the last 
200 years. We sampled mirids across both arid and agricul- 
tural regions in Australia, covering most of the geographic 
distribution of this endemic species and including the 
major host plants. We genotyped microsatellites from sam- 
ples spanning two seasons and sequenced a mitochondrial 
COI fragment from green mirids collected over 24 years. 
Given the ecology of green mirids outlined below, and the 
challenges posed by both agricultural and arid environ- 
ments, we structured our analyses according to the follow- 
ing three questions: (1) Do seasonal fluctuations in 
population size in both arid and agricultural regions result 
in genetic signatures of bottlenecks and drift? (2) Does 
long-distance dispersal occur between arid and agricultural 
populations? (3) Is genetic differentiation associated with 
the use of multiple host plants in the (arid) native range? 
We found that genetic patterns in C. dilutus are temporally 
dynamic, consistent with spatial and temporal heterogene- 
ity in its arid range. Long-distance dispersal between arid 
and agricultural populations is evident from the data, and 
host-associated differentiation was found between the pri- 
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mary host plants and alternative hosts in arid regions. 
Together, these results highlight the importance of consid- 
ering ecological and evolutionary processes across the dis- 
tribution of an organism. 

Materials and Methods 
Study system 

Creontiades dilutus is a mirid bug that is endemic to Aus- 
tralia and has become a major pest of several agricultural 
crops (Malipatil and Cassis 1997; McColl et al. 2011). 
This species has been associated with numerous host 
plant species and, prior to the advent of agriculture in 
Australia, was probably restricted to the relatively open 
interior. Here, the temporal variability of rainfall events is 
higher than in most other globally comparable desert sys- 
tems with similar mean annual rainfall (Morton et al. 
2011). Variability is compounded by years of drought 
(Nicholls 1991; Letnic and Dickman 2006). The availabil- 
ity of herbaceous plants is consequently ephemeral and 
often spatially patchy. During the summer months of 
November to February, C. dilutus persists in this region 
in low numbers despite temperatures in excess of 45°C 
and the availability of few host plants (JPH pers. obs., 
January 2007). Host plants persist longer in winter (June- 
August), but only if sufficient rain falls. With the excep- 
tion of "flood years," when host plants may be unusually 
widespread and persistent, suitable hosts generally require 
two or more successive rain events to thrive. These rain 
events usually occur only locally and typically the hosts 
are spatially patchy, with large areas of barren land 
between. Inland temperatures are close to optimal for 
mirid development during winter, allowing a generation 
time of around 25 days (Khan et al. 2009). Rapid popula- 
tion expansion is thus possible and large numbers of C. 
dilutus can be found where conditions are suitable. 

Creontiades dilutus presumably relies primarily on dis- 
persal to cope with adverse conditions by locating suitable 
patches of host plants when local conditions become 
unfavorable, because diapause (based on current evi- 
dence) is a facultative winter reproductive phenomenon 
(Miles 1995). Although C. dilutus persists throughout the 
year in the arid interior and sub-coastal agricultural 
regions of Australia, its abundance is seasonally inverse 
between the two, but the possibility of long-distance dis- 
persal remains untested. In inland Australia, C. dilutus 
relies on host plants that are both spatially and tempo- 
rally highly variable, even within a single season. Local 
population extinctions and founder effects might be 
expected, especially during dry years, as local resources 
die off and new patches are located. Conversely, when 
inland Australia experiences floods, the increased host 



abundance, together with the short generation time of 
this multivoltine insect, is likely to allow massive increases 
in population size across large areas. Pesticides applied in 
agricultural regions also have the potential to cause local- 
ized population contractions, and although green mirids 
are present on lucerne throughout the year in agricultural 
regions, their abundance during winter months is low 
(Miles 1995). 

Mitochondrial DNA sequencing and analysis 

Ten population samples were taken during 2006 and 2007 
from across Australia, and 146 of these individuals were 
sequenced for the COI fragment using the Folmer primers 
LCOI490 and HC02198 (Folmer et al. 1994) and a standard 
PCR protocol with an annealing temperature between 47°C 
and 50°C. We also obtained pinned specimens retained at 
The University of Queensland from previous research on 
this species. We were able to amplify the same fragment 
from 16 individuals collected in Gatton (Queensland) from 
lucerne in 1983 and 25 individuals collected in Biloela 
(Queensland) from lucerne in 1993. Details of the sample 
locations and Genbank accessions are provided in Table 1. 
DNA was extracted from the pinned specimens using 
DNeasy Blood and Tissue kit (Qiagen, Hilden, Germany) 
after being soaked in TE buffer overnight. The PCR proto- 
col was the same as for the ethanol preserved specimens. 
The COI fragments were sequenced bi-directionally at 
Macrogen (Korea) on an ABI3730, and then aligned, edi- 
ted, and trimmed using Codon Code Aligner v4.0 (Codon- 
Code Corporation, Dedham, MA). 

Haplotype networks were constructed using the R 
package TempNet (Prost and Anderson 2012). The tem- 
poral haplotype network was restricted to sites in the 
eastern Queensland cropping region (48 samples from 
2006/2007) where the samples from 1983 and 1993 had 
been collected. Nucleotide and haplotype diversity were 
calculated in DnaSP v. 5 (Librado and Rozas 2009). 

Microsate Mites - sample collection and 
genotyping 

A total of 32 population samples were collected from 17 
different host plant species in inland Australia and sub- 
coastal eastern Australia between January 2007 and March 
2008 (Table 2). Individual insects were preserved in 96% 
ethanol. DNA was extracted using a modified salt precipi- 
tation protocol based on that of Miller et al. (1988). Nine 
microsatellites (mirsat-2F, mirsat-4B, mirsat-3E, mirsat- 
Al, mirsat-3H, mirsat-6B, mirsat-5C, mirsat-G8, and 
mirsat-7G) were PCR amplified and genotyped on a 
Megabace capillary electrophoresis system (Amersham 
Biosciences) as per (Andris et al. 2010). Microsatellite 
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Table 1. Sampling locations for the 187 individuals for which the COI fragment was sequenced, including the number of sequences analyzed (A/) 
and the Genbank accession numbers. 



Location 


Date 


Lat (S) 


Long (E) 


Host plant 


N 


Genbank Accessions 


Adelaide 


2/12/2006 


-34.8208 


138.8700 


Polygonum convolvulus 


8 


JX186015 to JX186022 


Balingup 


14/09/2007 


-33.7889 


115.9760 


Solanum nigrum 


8 


JX1 86023 to JX1 86030 


BarcLong 


1 6/08/2006 


-23.5332 


145.0765 


Cullen cinereum 


8 


JX186031 to JX186038 


Biloela 


10/01/2007 


-24.3739 


150.5130 


Gossypium hirsutum* 


10 


JX1 86039 to JX1 86048 


Biloela 


10/01/2007 


-24.3739 


150.5130 


Medicago satlva* 


8 


JX1 86049 to JX1 86056 


Emerald 


14/08/2006 


-23.4958 


148.1884 


Verbesina enceliodes 


8 


JX1 86057 to JX1 86064 


Emerald 


1 5/08/2006 


-23.5722 


148.1001 


Verbesina enceliodes 


4 


JX1 86065 to JX1 86068 


Emerald 


1 5/08/2006 


-23.4663 


148.0918 


Vicia satlva 


8 


JX1 86069 to JX1 86076 


Kununurra 


28/08/2006 


-15.6459 


128.6969 


Gossypium hirsutum* 


5 


JX1 86077 to JX1 86081 


Longreach 


1 7/08/2006 


-23.4177 


144.2274 


Cullen cinereum 


8 


JX1 86082 to JX1 86089 


Longreach 


1 7/08/2006 


-23.4038 


144.2212 


Cullen cinereum 


8 


JX1 86090 to JX1 86097 


Longreach 


1 7/08/2006 


-23.4382 


144.2458 


Medicago polymorpha 


8 


JX1 86098 to JX186105 


Longreach 


1 7/08/2006 


-22.8941 


143.7867 


Swainsona galegifolia 


8 


JX186106 to JX186113 


Narrabri 


22/01/2007 


-30.2008 


149.5724 


Gossypium hirsutum* 


6 


JX186114to JX186119 


Narrabri 


22/01/2007 


-30.2008 


149.5724 


Medicago sativa* 


7 


JX186120 to JX186126 


Walget 


31/08/2006 


-29.9124 


146.9179 


Rapistrum rugosum 


8 


JX186127 to JX186134 


WintJun 


1 8/08/2006 


-22.4120 


143.0585 


Cullen cinereum 


7 


JX186135 to JX186141 


Wintlun 


1 9/08/2006 


-23.7810 


142.4658 


Cullen cinereum 


4 


JX186142 to JX186145 


WintJun 


1 9/08/2006 


-23.7338 


142.4287 


Senna artemisioides 


5 


JX186146 to JX186150 


Biloela 


20/09/1 993 


-24.3739 


150.5130 


Medicago sativa* 


12 


JX186151 to JX186162 


Biloela 


5/05/1993 


-24.3739 


150.5130 


Medicago sativa* 


13 


JX186163 to JX186175 


Gatton 


16/01/1983 


-27.5876 


152.3618 


Medicago sativa* 


16 


JX186176 to JX186191 


Byee 


14/03/2006 


-26.2566 


151.8539 


Cajanus cajan 


10 


EF016724 to EF016733 



'indicate samples from agricultural host plants. 



peaks were confirmed and binned manually. In total, 768 
specimens were genotyped; the DNA extractions of indi- 
viduals that failed to amplify at more than six loci were 
assumed to be low quality and were discarded, leaving 
665 genotyped individuals (Table 2), subsequent quantifi- 
cation of the samples with poor amplification revealed 
that the DNA was lower quality. 

HWE, genetic diversity, genetic 
differentiation, and tests for recent 
bottlenecks 

We estimated null allele frequency using the expectation 
maximization algorithm of Dempster et al. (1977) imple- 
mented in FreeNA (Chapuis and Estoup 2007) with 
10,000 bootstrap resamplings. Deviations from Hardy- 
Weinberg Equilibrium (HWE) were calculated using the 
exact probability test (Guo and Thompson 1992) imple- 
mented in Genepop (Rousset 2008) and a sequential Bon- 
ferroni correction was applied per locus to account for 
multiple tests. Locus mirsat-3E showed deviations from 
HWE in many samples and was consequently shown to 
have relatively high frequencies of null alleles (Table 3) 
and was discarded. The total number of alleles per locus, 
average number of alleles per locus, and (Nei 1987) unbi- 
ased gene diversity (per locus and sample) were calculated 



using FSTAT (Goudet 2001). Expected (He), observed 
(Ho), and unbiased expected (UHe) heterozygosities were 
computed using Genalex (Peakall and Smouse 2006). 
Exact tests for linkage disequilibria were carried out in 
Genepop (Rousset 2008). 

The proportion of genetic variance that can be attrib- 
uted to within-population comparisons and between-pop- 
ulation comparisons was estimated using an analysis of 
molecular variance (AMOVA) in Genalex (Peakall and 
Smouse 2006). Unbiased pairwise and locus specific F sr 's 
(Weir 1996) were computed with and without the algo- 
rithm for the exclusion of null alleles (ENA) implemented 
in FreeNA (Chapuis and Estoup 2007). Pairwise exact 
tests of genotypic differentiation were computed using 
Genepop (Rousset 2008), as this estimator is more appro- 
priate in situations where gene frequencies may deviate 
from HWE expectations, and a sequential Bonferroni 
adjustment was performed to account for multiple popu- 
lation comparisons. 

To test for signatures of recent demographic bottle- 
necks in the microsatellite data, the Wilcoxon test for het- 
erozygote excess (under the two-phase mutation model) 
and the allele frequency mode shift analysis were per- 
formed using the program BOTTLENECK (Piry et al. 
1999) for all 32 populations, and a sequential Bonferroni 
adjustment applied. 
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Table 2. Population codes, number of individuals genotyped, collection details, and host plant species for population samples used in this study. 



Code 


N 


Location 


Date 


Lat (S) 


Long (E) 


Host Plant 


Family 


BIL-GH* 


19 


Biloeala 


9/01/2007 


-24.3739 


150.5130 


Gossypium hirsutum 


Malvaceae 


BIL-MS1 * 


16 


Biloeala 


9/01/2007 


-24.3739 


150.5130 


Medicago sativa 


Fabaceae 


BIL-MS2* 


15 


Biloeala 


28/07/2007 


-24.3770 


150.5231 


Medicago sativa 


Fabaceae 


EMR-VE1* 


24 


Emerald 


10/01/2007 


-23.5723 


148.1002 


Verbesina encelioides 


Asteraceae 


EMR-MS* 


29 


Emerald 


10/01/2007 


-23.5772 


148.2084 


Medicago sativa 


Fabaceae 


EMR-GH* 


26 


Emerald 


10/01/2007 


-23.4663 


148.0916 


Gossypium hirsutum 


Malvaceae 


EMR-CA* 


21 


Emerald 


29/07/2007 


-23.4657 


148.0923 


Cicer arietinum 


Fabaceae 


EMR-VE2* 


26 


Emerald 


29/07/2007 


-23.5558 


148.1139 


Verbesina encelioides 


Asteraceae 


BIR-BS 


12 


Birdsvi lie 


3/08/2007 


-26.6707 


139.0784 


Epaltes cunninghamii 


Asteraceae 


SIM-BP1 


19 


Simpson 


4/08/2007 


-26.5743 


137.2792 


Blennodia pterosperma 


Brassicaceae 


SIM-GC 


13 


Simpson 


4/08/2007 


-26.5743 


137.2792 


Goodenia cycloptera 


Goodeniaceae 


SIM-SG1 


13 


Simpson 


4/08/2007 


-26.5742 


137.2745 


Senecio gregorii 


Asteraceae 


SIM-BP2 


9 


Simpson 


5/08/2007 


-26.3232 


137.0422 


Blennodia pterosperma 


Brassicaceae 


SIM-SG2 


1 1 


Simpson 


5/08/2007 


-26.3232 


137.0422 


Senecio gregorii 


Asteraceae 


SIM-CE 


23 


Simpson 


6/08/2007 


-25.9011 


138.8192 


Crotalaria eremaea 


Fabaceae 


SIM-CA 


29 


Simpson 


7/08/2007 


-25.9011 


138.8192 


Cullen australasicum 


Fabaceae 


EYR-CA 


28 


Eyre Creek 


7/08/2007 


-25.9002 


138.8540 


Cullen australasicum 


Fabaceae 


MIL-TS 


29 


Milparinka 


1 0/08/2007 


-29.5757 


141.9153 


Trigonella suavissima 


Fabaceae 


MIL-SG 


26 


Milparinka 


1 0/08/2007 


-29.5757 


141.9153 


Swainsona galegifolia 


Fabaceae 


MIL-SI 


29 


Milparinka 


10/08/2007 


-29.5757 


141.9153 


Sisymbrium irio 


Brassicaceae 


MIL-CP 


17 


Milparinka 


1 0/08/2007 


-29.5757 


141.9153 


Cullen pallidum 


Fabaceae 


TIL-EC 


26 


Tilpa 


1 1/08/2007 


-30.9365 


144.4160 


Erodium cygnorum 


Geraniaceae 


BOU-EC 


31 


Bourke 


1 1/08/2007 


-30.1743 


145.8144 


Erodium cygnorum 


Geraniaceae 


WAL-MP 


26 


Walget 


1 1/03/2008 


-30.0113 


148.0635 


Malva parviflora 


Malvaceae 


BRE-MP* 


19 


Brewarrina 


1 1/03/2008 


-29.9514 


146.7633 


Malva parviflora 


Malvaceae 


BRE-MS1 * 


25 


Brewarrina 


1 1/03/2008 


-29.9621 


146.8505 


Medicago sativa 


Fabaceae 


BRE-EC* 


26 


Brewarrina 


1 2/08/2007 


-29.9520 


146.7766 


Erodium cygnorum 


Geraniaceae 


BRE-PC* 


10 


Brewarrina 


1 2/08/2007 


-29.9520 


146.7766 


Phlegmatospermum cochlearinum 


Brassicaceae 


BRE-MS2* 


12 


Brewarrina 


1 2/08/2007 


-29.9615 


146.3504 


Medicago sativa 


Fabaceae 


NAR-MP* 


9 


Narrabri 


1 3/08/2007 


-30.3256 


149.7922 


Malva parviflora 


Malvaceae 


NAR-GH* 


24 


Narrabri 


20/01/2007 


-30.2008 


149.5724 


Gossypium hirsutum 


Malvaceae 


NAR-MS* 


23 


Narrabri 


20/01/2007 


-30.2008 


149.5724 


Medicago sativa 


Fabaceae 


♦indicate sampl 


es from agricultural sites. 












Table 3. Locus specific 


details for microsatellites used 


in this study; XWa, 


total number of alleles, Na, average number of al 


leles per population 


sampled, Ho, observed, and He, expected heterozygosities, HWD (number of population 


samples deviating from Hardy-Weinberg Equilibrium), 


null allele frequencies, and locus specific 


global F S t without and with the exclusion of null alleles. 




Locus 




Na 


Ho 


He 


HWD 


Null alleles gF ST Null 


gF ST No Null 


mirsat-2F 


11 


4.88 


0.52 


0.56 


2 


0.051 0.17 


0.17 


mirsat-4B 


10 


3.61 


0.31 


0.46 


2 


0.080 0.27 


0.25 


mirsat-3E 


16 


5.97 


0.25 


0.56 


13 


0.199 0.21 


0.18 


mirsat-A1 


9 


2.52 


0.19 


0.26 


5 


0.074 0.19 


0.21 


mirsat-3H 


21 


5.48 


0.41 


0.42 


1 


0.048 0.09 


0.08 


mirsat-6B 


6 


1.76 


0.03 


0.09 


0 


0.067 0.03 


0.09 


mirsat-BC 


20 


7.94 


0.82 


0.77 


0 


0.034 0.04 


0.04 


mirsat-G8 


13 


4.67 


0.30 


0.54 


6 


0.174 0.03 


0.03 


mirsat-7G 


15 


5.91 


0.69 


0.68 


2 


0.038 0.06 


0.06 



Spatiotemporal patterns of genetic 
differentiation and tests for migration 

We tested for the presence of isolation by distance (IBD) 
to explore gene flow in relation to the temporal and geo- 



graphic aspects of the sampling strategy. Initially, this 
analysis was restricted to the sampling period of July to 
August 2007 (when samples were widespread geographi- 
cally and collected over a short period from both agricul- 
tural hosts and native inland hosts). Subsequently, all 
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samples were analyzed to assess the temporal stability of 
the August 2007 pattern, these additional samples repre- 
sented the same agricultural crops sampled in January 
2007. The presence of an IBD effect was investigated by 
regressing ENA corrected genetic distance CFst/(1 — -Fst)) 
against geographic distance (Rousset 1997). A Mantel test 
of matrix correspondence was used to test for significance 
using the Isolation By Distance Web Service (IBDWS) 
3.15 (Jensen et al. 2005). 

Patterns of genetic differentiation and admixture, 
which may be obscured by statistics that assume the cor- 
rect a priori identification of populations, were clarified 
with the individual-based Bayesian clustering algorithm 
implemented in the program STRUCTURE (Pritchard 
et al. 2000). Low levels of null alleles are unlikely to affect 
the overall outcome of assignment testing such as the one 
implemented in the STRUCTURE algorithm (Carlsson 
2008). The "admixture" model was used as the most 
appropriate for a species in which dispersal is likely. Ini- 
tially, values of K from 1 to 14 were used, with a burn-in 
of 50,000 and a run length of 500,000, and each K value 
was replicated three times. These results were exported to 
the program STRUCTURE HARVESTER (Earl and Von- 
Holdt 20 1 1 ) and the most likely value of K for the data 
set was inferred using the AK method of Evanno et al. 
(2005). The data were then analyzed using this value for 
K with a burn-in of 100,000 and 1,000,000 subsequent 
iterations; this was replicated 10 times. The results were 
permuted with CLUMPP (Jakobsson and Rosenberg 
2007) and the mean of the permuted results plotted using 
DISTRUCT (Rosenberg 2004). 

Recent migration between arid inland sites and the 
eastern cropping regions was tested with the Bayesian 
assignment-based algorithm implemented in BAYESASS, 
which estimates rates of recent migration (within a few 
generations) from the proportion of individuals in each 
population sample that are assigned to other populations 
with high probability (Wilson and Rannala 2003). This 
algorithm does not assume that each designated popula- 
tion is in Hardy-Weinberg equilibrium, and produces 
reasonable estimates of actual migration in an experimen- 
tal setting (Mardulyn et al. 2008). Furthermore, simula- 
tions indicate that the method provides fairly accurate 
estimates of m if genetic differentiation is not too low 
(fsT > 0.05) (Faubet et al. 2007), this is the case for all 
of the between-location comparisons used in our analyses. 
The results of the IBD test indicated that genetic differen- 
tiation is unstable across seasons, so we restricted our 
analysis of gene flow to estimators of recent migration 
rather than coalescent approaches to estimating long-term 
averages of migration such as MIGRATE (Beerli and Fel- 
senstein 2001). Population samples were grouped by loca- 
tion and then split by the time of collection, because the 
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results of the STRUCTURE analysis indicated significant 
temporal shifts in cluster assignment within sites. We ran 
the analyses for 3,000,000 iterations with a burn-in of 
1,000,000 using the default delta value of 0.15 for allele 
frequency, migration rate, and inbreeding; the total log 
likelihood was plotted to assess convergence within runs. 
The computation was then performed 10 times with dif- 
ferent starting seeds to assess convergence across runs by 
comparing the posterior probability densities of allele fre- 
quencies. The results of the 10 runs were converted to 
tabular format using a custom Perl script (available on 
request from the corresponding author) for comparison. 
The number of times each outcome was achieved over 
the 10 runs was recorded, and the mean migration rates 
were calculated for each of these outcomes. Migration 
rates with lower 95% confidence intervals below m = 0.02 
were not considered significant and were also omitted. 
We used the lower 95% CI to assess the significance of 
migration rates because experimental tests in Caenorhabditis 
retnanei indicate that actual migration rates tend to be 
lower than the inferred rates, but within the 95% CI 
(Mardulyn et al. 2008). 

Host plant-associated differentiation 

Hierarchical AMOVA was performed across all 32 popu- 
lations, with the higher order defined as host plant in 
Genalex (Peakall and Smouse 2006). The host plants in 
central Australia, however, are completely different species 
to those used by green mirids in eastern cropping regions. 
To test the hypothesis that differentiation might be asso- 
ciated with host plant usage, we tested for genetic differ- 
entiation across host plants with respect to two inland 
localities within which multiple host plants had been 
sampled (namely Simpson Desert and Milparinka). The 
STRUCTURE algorithm was run using the admixture 
model with a burn-in of 100,000 and 500,000 subsequent 
iterations, with K = 2 "population" clusters. The most 
likely value of K was determined by the delta K method 
(Evanno et al. 2005) implemented in STRUCTURE HAR- 
VESTER (Earl and VonHoldt 2011). 

Results 
Mitochondrial 

Genetic diversity was low for the mitochondrial COI 
sequences, with one haplotype dominating each of the 
three temporal samples. The 2006/2007 samples (n = 146) 
had a nucleotide diversity (71) of 0.00055, and a haplotype 
diversity of 0.278. Tajima's D was -2.26 (P < 0.01) indi- 
cating an excess of low frequency polymorphisms likely 
due to population expansion. This pattern was similar in 
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1983 (n = 16), where nucleotide diversity was 0.00021, 
haplotype diversity 0.125, and Tajima's D —1.16, but not 
significant (P > 0.10), and in 1993 (n = 25), n = 0.00026, 
haplotype diversity 0.153, and Tajima's D —0.69, also not 
significant (P > 0.10). 

In the eastern cropping regions of Queensland, where a 
comparison of haplotype frequencies could be made 
across three temporal samples spanning 24 years, there 
has been a shift in the dominant haplotype between 1993 
and the more recent samples (Fig. 1, bottom). The haplo- 
type that was prevalent across the whole of Australia in 
2006/2007 (Fig. 1 middle) was present at much lower fre- 
quency in the earlier samples (Gatton, 1983 and Biloela 
1993) where the same haplotype was predominant despite 
a distance of 400 km between these two sites. The Biloela 
samples (1993 and 2006/2007) were all collected from the 
Queensland Government Biloela Research station 
(Table 1). 

HWE, genetic diversity, genetic 
differentiation, and tests for recent 
bottlenecks 

A total of 105 alleles weresscored across all loci and all 
populations, once the null allele prone locus mirsat-3E 
had been removed (Table 3). Unbiased gene diversity for 
each population (Nei 1987), when averaged across loci, 
ranged between 0.32 and 0.79 (mean 0.51) and was not 
significantly different between samples from agriculture 
and those taken inland (Fig. 2, two-tailed permutation 
test, P = 0.275). Four of 32 population samples, three 
from inland and one from agriculture, showed genetic 
signatures indicating a recent bottleneck in the allele 
mode shift analysis (BIR-BS, SIM-SG1, SIM-SG2, BIL- 
GH, Fig. 2), although only two of these showed a signifi- 
cant heterozygote excess in the Wilcoxon test (BIR-BS, 
P = 0.0117 and BIL-GH, P = 0.0078). In addition, three 
of the four populations showed indications of admixture, 
and neither of the Wilcoxon tests was significant (with an 
alpha probability of P > 0.05 after Bonferroni correction). 

Deviations from HWE were inferred in all loci for some 
populations, and the presence of null alleles was also 
inferred (Table 3). HWE deviations might, however, be 
expected in recently admixed populations due to the Wahl- 
und effect. We took three approaches to assess and mini- 
mize the effects of null alleles: 1) estimation of _F SX values 
using null allele corrected and non-corrected data; 2) 
removal of the two loci that had the greatest effect on HWE 
(mirsatlAl and mirsat2G8), and then comparing results 
across the 6 locus and 8 locus data sets; and 3) selection of 
analyses that are more robust to low frequencies of null 
alleles and small deviations from HWE (see methods for 
details). Evaluating the effects of null allele/HWE deviations 
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Figure 1. Map (top) showing sampling locations for the COI 
sequences obtained for this study (location colors maintained 
throughout the figure). Haplotype network (middle) showing all 
sequences from 2006/2007, and temporal haplotype network 
(bottom) showing the change in haplotype frequencies in Eastern 
Queensland between 1983 and 2006/2007. The area of each 
haplotype (circles and ellipses) represents the number of individuals 
having that haplotype (numbers inside haplotypes), empty ellipses 
(bottom figure) show haplotypes present in the 2006/2007 samples, 
but not in 1993 or 1983. 

using these three methods revealed that the low frequencies 
of null alleles inferred in some population samples for some 
loci did not dramatically affect the overall signal in the data, 
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Figure 2. Nei's unbiased gene diversity averaged across loci for all 
populations. Gray boxes represent samples collected from agricultural 
crops; black boxes represent samples collected from non-crop hosts. 
Asterisks indicate population samples for which there is some 
evidence of a recent bottleneck (see Results for details). 



and all results shown are for the 8-locus data set. Tests for 
genotypic linkage disequilibria returned no significant asso- 
ciations between pairs of loci for any of the 32 population 
samples after sequential Bonferroni correction for multiple 
comparisons. 

The AMOVA apportioned 19% of all molecular vari- 
ance to among-population comparisons (<5 PT = 0.188, 
P = 0.001). The global F ST estimates were similar with or 
without the elimination of null alleles, with the uncor- 
rected data returning only a slightly higher estimate 
(F ST = 0.122 using the ENA algorithm and 0.128 without 
ENA correction). Pairwise Fst's ranged from 0.0019 to 
0.329 (mean = 0.112), with 374 of 528 comparisons of 
genotypic differentiation being significant after sequential 
Bonferroni correction (Table SI). 

Spatiotemporal patterns of genetic 
differentiation and tests for migration 

The Mantel test of correspondence between geographic 
distance and genetic differentiation revealed a significant 
isolation by distance pattern when the analysis was 
restricted to the broad-scale geographic sampling of July 
to August 2007 (Fig. 3, r = 0.2897, P = 0.0099). In con- 
trast, when all sampling events were included in the anal- 
ysis (January 2007-March 2008), no isolation by distance 
effect was evident (r = 0.0076, P = 0.4465). The inclusion 
of these additional samples represented the same host 
plants that were sampled in agricultural regions during 
the July-August 2007 collections. 

The AK method (Evanno et al. 2005) indicated that 
K = 3 was the most likely number of genetic clusters for 
this data set. The combined and permuted results of the 
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Figure 3. Results of Mantel test for isolation by distance. Above: 
Samples collected during August 2007 (r = 0.2897, P = 0.0099). 
Below: All samples collected (January 2007-March 2008) (r = 0.0076, 
P = 0.4465). Solid line shows the matrix correlations, and dashed lines 
show 95% confidence intervals. 



subsequent 10 runs of the STRUCTURE algorithm with K 
set at 3 are shown in relation to the geographic origin of 
the population samples (Fig. 4). Individuals from Milpar- 
inka, Tilpa, and Bourke are mostly assigned to one cluster 
with high posterior probabilities, and these populations 
yielded the highest pair-wise F ST values when compared 
with the other sites (F sr = 0.047-0.307; mean = 0.132; 
148 of 156 tests of genotypic differentiation significant). 
Admixture was evident in several populations, as several 
individuals had mixed posterior probabilities of assign- 
ment to each of the other two clusters. The proportion of 
admixed individuals and their cluster assignment shifted 
between January 2007 and July 2007 in both Biloela and 
Emerald, and between January 2007 and August 2007 in 
Narrabri (Fig. 4.) 

We detected significant levels of recent migration, at 
the full geographic extent of sampling, using the BAYE- 
SASS algorithm. Variability was detected across runs of 
the algorithm, but this was characterized as a reversal of 
the inferred direction of migration between sites rather 
than changes in the sites between which dispersal was 



814 



) 2013 The Authors. Published by Blackwell Publishing Ltd. 



J. P. Hereward et al. 



Gene Flow in the Arid Adapted Green Mirid 




inferred. Figure 5 is a graphical representation of the 
migration rates and the frequency that each migration 
outcome was reached over 10 runs of the algorithm using 
different starting seed (See Table S2 for the full results, 
including 95% CI's). Significant migration was inferred 
between the Simpson Desert sites, in the arid interior, 
and the sub-coastal agricultural areas in Queensland 
(Biloela and Emerald), for example, Simpson (August 
2007) to Biloela (January 2007), m = 0.14, lower 95% 
CI = 0.08, upper 95% CI = 0.21, 6/10 runs. Although the 
direction of migration was most often toward agricultural 
regions, the direction was not always consistent across 
runs and strong inference cannot be made as to the direc- 
tion of dispersal from this result. Migration between Mil- 
parinka and other sites was not inferred from the genetic 
data, which is consistent with the outcomes of the 
STRUCTURE analysis and pairwise Fst's. 

Host plant-associated differentiation 

The Hierarchical AMOVA indicated an effect of host 
plant on molecular variance (^ RT = 0.078, P = 0.001); 
however, host plant species were not sampled consistently 
across the whole of the sampling area (because each has a 
restricted distribution relative to the scale of the study). 
Genetic differentiation in relation to host plant species was 




171 = 0.02-0.075 
m= 0.075 -0.125 
m= 0.125 -0.175 
m = 0.175-0.225 



Figure 5. Graphical representation of migration rates inferred using 
the Bayesian assignment algorithm in BayesAss. The size of the 
arrows indicates the migration rate (m), whereas the shade of the 
arrows indicates the number of times this outcome was reached over 
10 runs with varying starting seeds. Outcomes that were reached less 
than four times, and migration rates lower than 0.02 are not shown. 
(See text and Table S2, supporting information, for 95% CI's). 
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Figure 6. STRUCTURE analysis restricted to sites where multiple host plants were present and sampled. Each bar represents the posterior 
probability that the individual belongs to one of two clusters. The host plant species from which the particular samples were collected is listed 
below each population. 



therefore evaluated across two sites where several species 
could be sampled at each. At both of these, some degree 
of genetic differentiation was associated with plants in the 
genus Cullen relative to all other host plants sampled, 
although this was more pronounced at the Simpson Desert 
sites than at Milparinka in western New South Wales 
(Fig. 6) (but note that at Milparinka, C. dilutus was sam- 
pled from Cu. pallidum and at the Simpson Desert/Eyre 
creek sites from Cu. australaskum). The Evanno et al. 
(2005) delta K method indicated that K = 2 was the most 
likely number of clusters for the Simpson site, but there was 
no peak in delta K for the Milparinka site. Further examina- 
tion of the genotyping data at these two sites revealed that 
the genetic differentiation indicated by the STRUCTURE 
analysis appears to stem from the higher occurrence of rare 
alleles on Cullen host plants than on alternative hosts. 

Discussion 

Our combined analyses of genetic patterns in green mirid 
populations indicate that they are likely shaped by 
changes in population size (Fig. 1), admixture (Fig. 4), 
and recent migration (Fig. 3, Fig. 6). Our results are con- 
sistent with the expectations of complex spatiotemporal 
dynamics (Whitlock 1992; Wegmann et al. 2006). These 
dynamics are likely to arise from the extreme spatial and 
temporal environmental heterogeneity that typify dry 
years in this bug's native range in arid Australia. Hetero- 
zygote excesses and allele frequency mode shifts indicate 
that several populations may have passed through recent 
localized bottlenecks, although alternate processes can 



cause heterozygote excess. However, the prevalence of a 
single COI haplotype across the whole of the continent 
(from 2006/2007 samples), negative Tajima D values, low 
haplotype diversity in the older samples, and the shift in 
the most prevalent haplotype between 1993 and 2006/ 
2007 (Fig. 1) are consistent with recurrent reductions in 
population size over a longer period. These past reduc- 
tions in population size may reflect alternate periods of 
drought (when resources are ephemeral and spatially pat- 
chy) and floods, which cause widespread environmental 
homogeneity in terms of host plant availability. The pres- 
ent broad geographic prevalence of the one COI haplo- 
type, inference of recent migration from the microsatellite 
data (Fig. 5), and admixture in agricultural populations 
(Fig. 4) demonstrate that populations in arid and agricul- 
tural regions are connected by gene flow. Although weak 
genetic differentiation was detected locally (within arid 
regions) across their primary host plants (in the genus 
Cullen) and alternative plant species, it was mostly the 
presence of rare alleles that was responsible for this pat- 
tern. We suggest that rare alleles on the primary hosts 
(Cu. cinereum and Cu. australasicum) could be a conse- 
quence of higher abundance of C. dilutus relative to their 
abundance on alternative host species. These points are 
expanded and justified below. 

Genetic diversity and founder effects 

We assessed genetic diversity and possible founder effects 
across inland and agricultural sites to determine whether 
ephemeral host availability (inland) or pesticide use (in 
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agriculture) affected the temporal stability of patterns of 
genetic differentiation between mirid populations. We 
found no significant difference in microsatellite diversity 
between C. dilutus collected from crop hosts and those 
collected from non-crop hosts (Fig. 2). Genetic signals of 
recent bottleneck events (heterozygote excess) were pres- 
ent in the microsatellite data in three inland and one 
agricultural population of C. dilutus (Fig. 2). None of the 
tests, however, were significant after Bonferroni correc- 
tion, so we can only tentatively infer localized contrac- 
tions. Lucerne, the primary crop host of C. dilutus, is 
often grown without pesticides in Australia, and insecti- 
cide-induced bottlenecks are less likely to be driving pat- 
terns of genetic differentiation than the arid dynamics of 
this mirid species. 

Genetic differentiation among C. dilutus populations 
was higher than generally reported for pest insect species 
surveyed within an agricultural context (Endersby et al. 
2006, 2007; Kim et al. 2009; Torres and Azeredo-Espin 
2009). The highest F ST values were attributed to pair-wise 
comparisons among three arid inland sites in New South 
Wales (Milparinka, Tilpa, and Bourke), and other popula- 
tions. The Structure analysis also clearly differentiated 
populations at these three sites from others. Although no 
heterozygote excess was detected at these sites (indicative 
of a recent bottleneck), we suspect that a combination of 
spatial heterogeneity and founder effects could contribute 
to strong genetic drift (and therefore high genetic differ- 
entiation) given that our study was conducted during a 
dry period when patches of host plants were separated by 
large areas of barren land. Elevated Pst's are predicted 
(even when migration rates are high) under spatiotempo- 
rally dynamic population models if environmental hetero- 
geneity contributes to a large variance in local population 
size (Wegmann et al. 2006), a scenario consistent with 
the ecology of C. dilutus. 

The low nucleotide diversity in both the older and 
more recent samples (2006/2007 n = 0.00055, 
1983 = 0.00021, 1993 = 0.00026) and change in predomi- 
nant mitochondrial hap lo type over the last 13 years was 
striking. Mitochondrial DNA is expected to suffer a more 
extreme loss of alleles than nuclear markers during demo- 
graphic bottlenecks due to the uniparental inheritance of 
the plastid and the reduced effective population size of its 
genome (Wilson et al. 1985; Simon et al. 1994). For 
example, a local population founded by a single gravid 
female would have one mitochondrial haplotype, but 
potentially four microsatellite alleles. Similar shallow "star 
shaped" genealogies and negative values of Tajima's D 
have been reported in agriculturally damaging insects 
with documented dispersal capacity (Albernaz et al. 
2012), and in the case of the widespread noctuid pest 
Helicoverpa armigera, this pattern even spans continents 



(Behere et al. 2007). Human assisted range expansion of 
pest insects through the provision of agricultural 
resources is the scenario that typically explains widespread 
haplotypes (Grapputo et al. 2005). In green mirids, how- 
ever, the dominant haplotype not only occurs across both 
agricultural and native arid regions but also has changed 
within the last 24 years, indicating that the alternation 
between dry and wet years in arid regions could be 
responsible, rather than the introduction of agricultural 
resources over the last 200 years. The low haplotype 
diversity in the older samples indicates that at least two 
such population contractions are likely to have occurred. 

Long-distance dispersal between arid and 
agricultural populations 

That one haplotype is now dominant across the 5000-km 
width of Australia indicates that dispersal in C. dilutus 
has been widespread. Geographic differentiation was 
higher in the microsatellite data set, but the geographic 
distribution was not stable over time. This is evident from 
the temporary nature of the isolation by distance effect 
(Fig. 3), the temporal shifts in cluster assignment in the 
STRUCTURE analysis at Biloela, Emerald, and Narrabri 
(Fig. 4), and by the change in most prevalent COI haplo- 
type between 1993 and 2006. Admixture across large geo- 
graphic distances most likely results from dispersal, and 
this is evidenced by the inference (B AYES ASS; Fig. 5) of 
significant migration rates across distances over 1500 km. 
The direction of inferred migration was not consistent 
across multiple runs of the algorithm, so conclusions 
regarding the directionality of dispersal remain tentative. 
Return migration from agricultural regions back to Cen- 
tral Australia by pest populations that derived originally 
from central desert areas (as postulated for Helicoverpa 
punctigera, which is also an Australian arid adapted spe- 
cies) is thought to be unlikely based on prevailing wind 
directions and because positive evidence of its existence 
has never been found (Downes et al. 2010). 

The sites between which migration was inferred were 
consistent across runs of the BayesAss algorithm, and are 
thus likely to represent sites between which recent migra- 
tion has occurred. Creontiades dilutus populations can 
expand rapidly, and abundance is seasonally inverse 
between inland and eastern regions. A migration event in 
late spring/early summer, when numbers are high in 
inland areas and low in cropping regions, might therefore 
result in a much higher inferred migration rate than the 
actual number of individuals migrating and establishing 
successfully. The regular seasonal influx of C. dilutus to 
cotton crops, which does not appear to be derived from 
local lucerne populations (Miles 1995), indicates that 
dispersal from inland populations may be a regular 
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occurrence, the microsatellite data support this hypothe- 
sis, but it does require further direct testing. 

Host plant-associated genetic 
differentiation in arid regions 

We found weak genetic differentiation between C. dilutus 
from Cu. australasicum and alternative hosts in the same 
geographic area (Simpson) in the structure analysis 
(Fig. 6), which may partly account for the significant role 
of host plants implicated by the hierarchical AMOVA. 
Plants in the genus Cullen maintain a significantly higher 
density of C. dilutus than other available hosts, indicating 
that plants in this genus are primary hosts for green mir- 
ids (Hereward and Walter 2012). However, analyses of 
gut contents using chloroplast intron markers revealed 
that a substantial proportion of C. dilutus individuals col- 
lected from the Cullen primary host plants had recently 
fed on other host plants (Hereward and Walter 2012). 
The use of multiple plant species by C. dilutus is perhaps 
best understood as a behavioral adaptation to survive in 
an arid environment where host plants are ephemeral and 
the primary host species not always available (Velasco 
and Walter 1993). The physiological and behavioral pro- 
cesses that underpin their multiple host use warrant 
investigation in association with the movement of indi- 
viduals. 

Previous quantified sampling showed that C. dilutus 
abundance was significantly higher on the Cullen hosts, 
Cu. cinereum and Cu. australasicum, than alternative hosts 
locally, but not Cu. pallidum (Hereward and Walter 
2012). We detect genetic differentiation in green mirids 
between Cu. australasicum and other hosts locally, but 
not for Cu. pallidum (Fig. 6), although this pattern 
requires further replication to ensure that it is not the 
result of fine scale spatial differentiation. More rare mi- 
crosatellite alleles were present in green mirid populations 
from Cu. australasicum than from alternative hosts. This 
may be a consequence of a much greater proportion of 
green mirids being attracted to these plants (perhaps from 
refuges provided by alternate hosts) and surviving. Fur- 
ther temporal samples at these inland sites would allow 
temporal estimates of effective population size across pri- 
mary and secondary host plants, and the relative absence 
(and perhaps even loss) of rare alleles on alternative hosts 
needs to be investigated directly if these patterns are to be 
understood mechanistically. 

Conclusions and implications 

Creontiades dilutus shows evidence of widespread dispersal 
in both the mitochondrial and microsatellite data sets 
examined here, despite this species using different plant 



resources (both locally and regionally), having seasonally 
inverse abundance between inland and agricultural 
regions, and presumably experiencing different selective 
pressures in these regions of Australia. The change in the 
most prevalent mitochondrial haplotype over 24 years is 
consistent with successive population contractions and 
expansions, likely in relation to fluctuations between dry 
periods and wet periods in the arid regions of Australia. 
Dispersal appears to be the major mechanism by which 
C. dilutus is able to survive on the ephemeral resources in 
this region, and the data provide no indication that large 
numbers of these bugs persist through dry periods by dia- 
pause. The spatiotemporal dynamics and changing gene 
frequencies outlined above contrast with the lack of 
genetic differentiation found in the same agricultural 
regions over several years for the highly dispersing H. ar- 
migera (Endersby et al. 2007) and the temporal stability 
of allele frequencies recoded for Queensland fruit fly (Yu 
et al. 2001). These dynamics also differ from the stepwise 
founder effects associated with insects that colonize new 
temporally stable habitat "islands" through human move- 
ment (Stone and Sunnucks 1993). Patterns of genetic dif- 
ferentiation and gene flow in green mirids seem to be 
driven instead by the spatial and temporal heterogeneity 
of their native hosts, but these same effects have spread 
to agricultural regions. This fits with Oliver (Oliver 
2006) 's hypothesis that the expansion of host resources is 
likely to increase gene flow in native insects. 

With such spatiotemporal dynamics, adaptation to 
novel host plants is unlikely. We nevertheless found weak 
host-associated differentiation between green mirids on 
their primary host plants and those on alternative hosts 
growing locally, despite establishing previously that these 
individuals will feed on alternative host species even when 
in the nearby vicinity of the primary host (Hereward and 
Walter 2012). Many herbivorous insects that use multiple 
hosts have been shown, by thorough quantitative sam- 
pling, to have a similar closer affinity to one host species 
than others that it may use (Milne and Walter 2000; 
Rajapakse et al. 2006; Manners and Walter 2009). Assess- 
ing gene flow and genetic diversity in many of these 
instances might further our understanding of multiple 
host use by herbivorous insects. 

Our results highlight the importance of assessing evolu- 
tionary and ecological processes across the distribution of 
an organism that uses both native and human-altered 
habitats simultaneously. If our analyses had been 
restricted to either agricultural areas or localized parts of 
the arid range of this species, our interpretations might 
be quite different. For example, broader geographic analy- 
ses of Rhagoletis pomonella, perhaps the most famous 
example of host-associated differentiation following the 
human introduction of novel hosts (cultivated apple) 
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(Bush 1993), to include native hosts in Mexico, indicates 
that the differences in host plant use had an allopatric 
rather than sympatric origin (Feder et al. 2003; Michel 
et al. 2007). In C. dilutus, we find that the spatiotemporal 
dynamics in its arid native range continue to drive genetic 
patterns across both arid and agricultural environments. 
The adaptations that allow it to persist despite the spatio- 
temporal heterogeneity of host resources in arid regions 
(migration and the use of alternative hosts) appear to 
have not only facilitated the colonization of new agricul- 
tural habitats but also maintain gene flow across large 
distances. 
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